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ABSTRACT 

Geometrically thick accretion flows may be present in black hole X-ray binaries observed in the 
low/hard state and in low-luminosity active galactic nuclei. Unlike in geometrically thin disks, the 
angular momentum axis in these sources is not expected to align with the black hole spin axis. 
We compute images from three-dimensional general relativistic magnetohydrodynamic simulations of 
misaligned (tilted) accretion flows using relativistic radiative transfer, and compare the estimated 
locations of the radiation edge with expectations from their aligned (untiltcd) counterparts. The 
radiation edge in the tilted simulations is independent of black hole spin for a tilt of 15°, in stark 
contrast to the results for untiltcd simulations, which agree with the monotonic dependence on spin 
expected from thin accretion disk theory. Synthetic emission line profiles from the tilted simulations 
depend strongly on the observer's azimuth, and exhibit unique features such as broad "blue wings." 
Coupled with precession, the azimuthal variation could generate time fluctuations in observed emission 
lines, which would be a clear "signature" of a tilted accretion flow. Finally, we evaluate the possibility 
that the observed low- and high-frequency quasi-periodic oscillations (QPOs) from black hole binaries 
could be produced by misaligned accretion flows. Although low-frequency QPOs from precessing, 
tilted disks remains a viable option, we find little evidence for significant power in our light curves in 
the frequency range of high-frequency QPOs. 

Subject headings: accretion, accretion disks — black hole physics — radiative transfer — relativity 

1. INTRODUCTION 

In standard t hin disk accretion theor y 
(jShakura fc Sunvaevi [l97l INovikov fc Thornel Il97l . 
the angular momentum axis of the accretion flow is 
assumed to be align e d with the black hole spin axis. 
iBardeen fe Pettersonl (|1975l) found that even if the 
initial angular momentum axis of the accretion flow 
is misaligned from the black hole spin axis, the inner 
part of the disk will still align on the viscous timescale. 
However, this so-called "viscous" regime only operates 
when H/R < a, where H/R is the scale height of the 
accretion disk, and a is the parameterized viscosity 
(|Papaloizou fc Lir] 119951 ) . This is applicable in active 
galactic nuclei (AGN) and the high/soft or thermal 
state of black hole X-ray binaries. On the other hand, 
advection-dominated accretion flows (ADAFs) are 
expected in the low/hard state of black hol e X-ray 
binaries (|Naravan fc Yil 119951: lEsin et al.l 119971 ) and in 
low-luminosity AGN. ADAFs are unable to cool through 
efficient radiation, and are geometrically thick. It is 
likely that the accretion flow in many of these sources is 
misaligned, or "tilted." 

Contempo rary general relativistic MHD simulations 
(GRM HD, IDe Villiers fc Hawlevl 12001 IGammie et all 
2003) currently provide the most physically realistic de- 
scription of the inner portion of accretion flows around 
spinning black holes. Radiation can be calculated 
from these simulations in post-processing by assuming 
that it is dynamically and thermodynamically negligible. 
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This method has been used to look for high frequency 
quasi-periodic oscillation s (HFQPOs) in simulated data 
(|Schnittman et al I2006D and to create radiative models 



of Sagittarius A* (INoble et al.ll2007t IMoscibrodzka et al.l 
2009: IDexter et al.H2009l 12010) . 



All of this work assumed alignment between the an- 
gular momentum ax is of the accretion flow a nd the 
black hole spin axis. iFraeile et al.l (I20071 . [2001) : iFraeild 
(2009) were the first to do GRMHD simulations of disks 
with a tilt between these two axes. These new simula- 
tions yielded a number of unexpected features. First, 
the main body of the disk remained tilted with re- 
spect to the symmetry plane of the black hole; thus 
there was no indication of a Bardeen-Petterson effect 
in the disk at large. The torque of the black hole 
instead principal ly caused a globa l precession of the 
main disk body (jFragile fc Anninosl 120051 : iFragile et al.1 
2007). The time-steady structure of the disk was also 
warped, with latitude-dependent radial epicyclic motion 
driven by pressure gr adients attributable to the warp 
(jFragile fc Blaesll2008| ). The tilted disks also truncated 
at a larger radius than expected for an untilted disk. 
In fact, based on dynamical measures, the inner edge of 
these tilt ed disks was found to be independent of black 
hole spin (|Fragilell2009D . in sharp contrast to the expecta- 
tion that accretion flows truncate at the margina l ly sta - 
ble orbit of the black hole. Finally, iHenisev et "all (2009) 
found evidence for trapped inertial waves in a simulation 
with a black spin a = 0.9, producing excess power at a 
frequency 118(M/10M Q )" 1 Hz. 
In this work we use relativistic ray tracing to produce 
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Figure 1. Sample images of the thermal emission model for the 
90h (left) and 915h (right) simulations at 60° inclination. The 
observed photon energy is Eo = lOkeV for a IOMq black hole, and 
each panel is 54M across. The color scale is linear, increasing from 
blue to red to yellow to white. 
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images and light curves of some of these numerically sim- 
ulated tilted and untilted black-hole accretion disks. Our 
goal in this paper is to discuss observable differences be- 
tween the two types of accretion flows, and to identify ob- 
servational signatures of tilted black hole accretion disks. 

2. METHODS 
2.1. Simulation Data 

Th e simulat i ons us ed here are from lFragile et all (2007, 
I2009T) : IFragile] (|2009D . The parameters are given in Table 
[TJ All of the simulations u sed the Cosmos++ GRMHD 
code (|Anninos et al.l 120051) . with an effective resolution 
of 128 3 for the spherical-polar grid (except near the 
poles where the grid was purposefully underresolved) and 
128 x 64 x 64 x 6 for the cubed-sphere grid. The simula- 
tions were initialized with an anal ytically solvable, time- 
stead y, axisymmetric gas torus (|De Villiers fe Hawlevi 
2003), threaded with a weak, purely poloidal magnetic 
field that follows the isodensity contours and has a min- 
imum Pgas I P-mag = 10 initially. The magnetorotational 
instability (MRI) arose naturally from the initial condi- 
tions, and the disk quickly became fully turbulent. The 
simulations were all evolved for ^8000M, or ~40 orbits 
at r = 10M in units with G = c = 1. Only data from 
the final 2/3 of the simulation are used in this analysis, 
once the disks are fully turbulent as measured by a peak 
in the accretion rate and in the mass inside of r = 10M. 
This is chosen to utilize as much of the simulation data 
as possible, and none of our results depend on which time 
interval in the simulation is used. 

These simulations all evolved an internal energy equa- 



tion, and injected entropy at shocks. Such a formulation 
does not conserve energy, and produces a more slender, 
cooler torus than conservative formulations which cap- 
ture t he heat from numerica l reconncction of magnetic 
fields (|Fragile fe Meierir2009() . The scale height spanned 
the range H/R ~ 0.05 — 0.1 in these simulations, with 
larger scale heights for higher spin simulations. 

2.2. Ray Tracing 

Relativistic radiative transfer is computed from simu- 
lation data via ray tracing. Starting from an observer's 
camera, rays are traced backwards in time assuming they 
are null geodesies (geometric opt ics approximati o n), us - 
ing the public code described in iDexter fc Agoll (|2009|) . 
In the region where rays intersect the accretion flow, the 
radiative transfer equation is solved alon g the geodesic 
(Brodcrick| I2006D in the form given in iFuerst fe Wul 
(2004), which then represents a pixel of the image. This 
procedure is repeated for many rays to produce an im- 
age, and at many time steps of the simulation to produce 
time-dependent images (movies). Light curves are com- 
puted by integrating over the individual images. Sam- 
ple images of two simulations are given in Figure [TJ 
Doppler beaming causes asymmetry in the intensity from 
approaching (left) and receding (right) fluid. Photons 
emitted from the far side of the accretion flow are de- 
flected toward the observer, causing it to appear above 
the black hole. The thick, central ring is due to gravita- 
tional lensing from material passing under the black hole, 
while the underresolved circular ring is caused by pho- 
tons that orbit the black hole one or more times before 
escaping. These ring features are in excellent agreement 
with the predictions made by iViergutzl (|1993l ). 

To calculate fluid properties at each point on a ray, 
the spacetime coordinates of the geodesic are trans- 
formed from Boyer-Lindquist to the Kerr-Schild coordi- 
nates used in the simulation. Since the accretion flow is 
dynamic, light travel time delays along the geodesic are 
taken into account. Data from the sixteen nearest zone 
centers (eight on the simulation grid over two time steps) 
are interpolated to each point on the geodesic. Between 
levels of resolution near the poles on the spherical-polar 
grid, data from the higher resolution layer are averaged to 
create synthetic lower resolution points, which are then 
interpolated. Very little emission originates in the un- 
derresolved regions of the simulation. 

The simulations provide mass density, pressure, veloc- 
ity and magnetic field in code units. These are con- 
verted into cgs un i ts follo wing the procedure descr ibed in 
iSchnittman et al.l ([20061 ) and IDexter et al.l (|2009| ). The 
length- and time-scales are set by the black hole mass, 
taken to be 1OM throughout. 

We consi der two emission models . The thin line emis- 
sivity from ISchnittman et all ()2006l ) is a toy model that 
traces the mass density in the accretion flow . The ther- 
mal emission model from ISchnittman et all (|2006f ) uses 
free- free emission and absorption coefficients, and is used 
as a model for the high/soft state. Although we do not 
expect tilted disks to accurately represent the high/soft 
state, this model may be appropriate for sources radiat- 
ing at an appreciable fraction of Eddington, where the 
infall time is shorter than the radiative diffusion time and 
the accretion flow becomes geometrically "slim." When 
taking the temperature from the ideal gas law rather 
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Figure 2. Comparison of relative intensities for all simulations 
using the thin line emissivity. The flux from grids of images over 
observer time, inclination and azimuth for each simulation have 
been averaged to create these curves. 



Figure 4. Comparison of relative intensities for all simulations 
using the thermal emissivity at Eo = 1 keV. The flux from grids of 
images over observer time, inclination and azimuth for each simu- 
lation have been averaged to create these curves. 
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Figure 3. Radiation edge as a function of spin for untilted (open) 
and tilted (solid) simulations for the thin line emissivity. The error 
bars show the one standard deviation time variability in the radi- 
ation edge, averaged over other parameters. The solid line is the 
marginally stable orbit. 

than the radiation-domi nated equation of state used in 
ISchnittman et al.l (|2006[) . this model may be qualita- 
tively appropriate for modeling the low/hard state in 
X-ray binaries or low-luminosity AGN. 

In £|3 - H we consider emission from inside of r = 15M, 
while in ^3.2! and I3.3I fluid inside of r = 25M is used for 
the ray tracing. For all results here, we take the tem- 
perature from the ideal gas law rather than assuming 
a radiation-dominated equation of state. All of our re- 
sults are qualitatively identical when using the radiation- 
dominated equation of state to calculate the tempera- 
ture. 

3. RESULTS 

3.1. Radiation Edge 

Inferring the inner edge of accretion flows is important 
for attempts to m easure spin from broad ir on lines (e.g. 
Wilms et al.l2001l) or con tinuum fitting (e.g. lShafee et al.l 
20061 iDavis et all 120061 ). Such measurements assume 
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Figure 5. Radiation edge as a function of spin for untilted 
(open) and tilted (solid) simulations for the thermal emissivity at 
Eo = IkeV. The error bars show the one standard deviation time 
variability in the radiation edge, averaged over other parameters. 
The solid line is the marginally stable orbit. 

that the disk has a sharp cutoff at the innermost sta- 
ble ci r cular orbit, wh ich depends on spin (Bardcc n et al.l 
I1972T) . IFragild (12009D used four diff erent dynamical mea- 
sures from lKrolik fc Hawlevl (|200l to compare the inner 
edges of simulated tilted and untilted accretion disks. 
Here we use the ray traced models to locate the "radi- 
ation edge," the radius inside of which the contribution 
to the total flux is negligible. 

For each emission model, images are calculated for all 
simulations over a grid of observer inclination, observer 
time, and observer azimuth (for the tilted simulations). 
We then compute images cutting out fluid inside of suc- 
cessive values of the radius r in . The radiation edge is 
functionally defined as the radius where the ratio of in- 
tensities, F(r in )/ F(0), drops below an arbitrary fraction 
/, chosen so that the untilted radiation edge agrees as 
well as possible with r ms , the marginally stable orbit. 

Figure [2] shows a plot of F(ri n )/F(0) as a function of 
ri n averaged over observer time, azimuth and inclination 
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for the thin line emissivity. From these curves we extract 
values of the radiation edge, r oc i go - Results are shown in 
Figure |3l where the error bars are computed from the 
standard deviation of r e( jge as a function of time, aver- 
aged over the other parameters. This result agr ees well 
with the dynamical measures from[Eagile (2009). While 
the radiation edge moves in towards the black hole with 
increasing spin for untiltcd simulations, there is no such 
trend in the tilted simulations. Instead, the radiation 
edge appears to be independent of spin. 

Figures |4] and [5] show the same plots for the thermal 
emission model with observed photon energy Eq = 1 
keV. The conclusions are identical with this emission 
model. The untilted simulations have radiation edges 
which agree quite well with r ms , while the tilted simula- 
tions show no correlation betwe en spin and r P ,d „ r .. Again, 
these results are consistent with [Fragile (2009), although 
we find no trend of increasing radiation edge with spin, 
as was found for a couple of the dynamical measures used 
in iFra gile (2009). Plots from other observed photon en- 
ergies are not shown; although the relative flux falls off 
much more quickly with increasing r; n at higher photon 
energies, the results for the radiation edge remain com- 
pletely unchanged. 

3.2. Emission Line Profiles 

Spectra from AGN and X-ray binaries typically include 
strong emission and absorption features. As the observed 
line shapes are sensitive to both the velocity of the emit- 
ting/absorbing fluid, and also to the local gravitational 
redshift, they can provid e information about the dynam - 
ics of the accretion flow ()Fabian et alj|1989t lLaorl(l99ll ). 

Untilted accretion flows have nearly Keplerian velocity 
distributions outside the marginally stable orbit, where 
the velocities smoothly transition to plunging. Simulated 
tilted accretion disks, on the other hand, show three ma- 
jor differences. The Keplerian velocity structure is now 
tilted. 

Secondly, the warped structure of the tilted disks leads 
to epicyclic motions with velocity magn itudes compara- 
ble to the local geodesic orbital velocity (jFragile &: Blaesl 
2008). Finally, the larger radiation edge values of the 
tilted disks identified in £13.11 means that the transition 
to plunging orbits occurs at larger radius than in untiltcd 
disks. 

These effects indicate that we should expect a number 
of differences in lin e profiles from tilted accretion flows 
(jFragile et all 12005). The maximum blueshift should be 
larger for tilted accretion disks, except for edge-on view- 
ing. For i < 90° — f3, where i is the observer's incli- 
nation angle and (3 is the initial tilt angle, both rela- 
tive to the black hole spin axis, the tilted accretion flow 
should mimic an untilted one with a larger inclination. 
In contrast, the red wing should be less pronounced in 
the tilted disks due to their larger truncation radii. On 
the redshifted side, tilted disks behave similar to lower 
spin, untilted disks. 

Producing a detailed reflection spectrum would require 
a significant number of assumptions to model the metal- 
licity, ionization levels, and incident X-ray flux through- 
out the accretion flow. For simplicity, we instead use 
toy model emissivities of the form j oc pr~ s , where p 
is the fluid mass density, j is the photon-energy inte- 
grated emissivity and s = 2,3. The two values corre- 



spond to assuming the emitted line flux is proportional 
to the incident flux from an irradiating source on the 
spin axis and to local d issipation of heat, respectively 
(e.g.. IFragile et aUl2005l ). This simple form allows us to 
focus on general features to be expected from emission 
lines from tilted black hole accretion disks. 

Figure |5] shows sample line profiles for an inclination of 
i = 60° for four observer azimuths from the 915h simula- 
tion. Only a single observer azimuth from the 90h sim- 
ulation is shown, since the time-averaged emission line 
is independent of observer azimuth for untilted simula- 
tions. In all cases, the lines consist of a strong peak near 
the rest energy of the line (g = E /E em — 1), a smaller 
peak at lower energy and a "red wing," whose extent and 
strength depends on the amount of emission arising very 
close to the black hole (small g). The location of the 
"blue" peak (large g) depends on the maximum veloc- 
ity along the line of sight in the accretion flow. For an 
untilted disk, this corresponds directly to the observer's 
inclination angle, since all fluid velocities are essentially 
in the equatorial plane. 

For the tilted model shown in Figure HI the loca- 
tion and strength of the blue peak changes significantly 
with observer azimuth. When the angular momentum 
axis of the accretion flow is in the plane of the sky 
(— tt/2 < 0o < — 7r/4, depending on the simulation time), 
its fluid velocities are maximally aligned with the ob- 
server's line of sight, leading to the largest blueshifts. 
This is the same condition as an untilted disk being 
viewed edge-on. For other orientations, the blue tail 
can extend to significantly higher photon energies in the 
tilted simulations because the largest effective inclina- 
tion is approximately i e s = i + ft- When the accretion 
flow is not edge-on, there will exist orientations where 
i e s > i, and the blue peak for a tilted simulation will oc- 
cur at higher energy than possible for untilted accretion 
flows. The red wing, on the other hand, remains largely 
unchanged with observer azimuth, since it is caused by 
gravitational redshifts rather than Doppler boosts. Since 
the radiation edge for the 915h simulation was found to 
occur at significantly larger radius than that of 90h, it is 
expected that the red wing should extend further in the 
90h simulation. The effect is subtle, but identifiable in 
Figure [H 

To quantify these trends, for all simulations we com- 
pute the extent of the line profile, as well as the strengths 
and locations of their red and blue peaks. Most clear are 
the results for the line extents, shown for i = 60° in 
Figures [7] and [5] As expected, the red wing extends to 
lower photon energies at higher spins for untilted simula- 
tions, while there is no similar trend for the tilted models. 
Also as expected, the blue wing extends to systematically 
higher photon energies in the tilted simulations because 
of the difference between i e g and i noted above and the 
epicyclic motion in the tilted simulations. 

Perhaps the most striking feature of the line profiles 
is the variation with observer azimuth seen in all tilted 
simulations. These changes in line shape between differ- 
ent observer azimuths are typically larger than the full 
range of changes seen between different spins for untiltcd 
simulations. This suggests that the most powerful means 
of recognizing a tilted accretion disk may be to measure 
changes in an emission line profile over time as the disk 
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Figure 6. Emission line profiles for simulations with a=0.9. The emissivity is j oc pr 3 and the observer inclination is 60° in all cases. 
The dotted lines show the la range, taken from the time variability. 




Figure 7. Minimum line energy vs. spin for all simulations. 
The tilted (untilted) simulations are denoted by solid (open) cir- 
cles, and four observer azimuths are plotted for the tilted simu- 
lations. The open diamonds are from a thin disk in the equa- 
torial plane with an emissi v ity j o c r~ 3 similar to that used in 
ISchnittman fc Bertschingerl |(2001 1; IDexter fc Keioj pOOl V The 
minimum line energy is defined as the lowest energy contained in 
the set of intensities comprising 99% of the total line intensity. The 
la errors are taken from the time variability. 
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Figure 8. As in Figure[7] but for the maximum line energy. 



precesses. 



3.3. Variability 



X-ray timing of black hole binaries has allowed 
the characterization of power spectra and the 
detection of transient QPO s (for a review, see 
iRemillard fc McClintockl I2006I ). High-frequency QPOs 
are seen in the steep power law state (SPL), while 
low-frequency QPOs have been observed in both the 
hard state and the SPL. The geometry of the accretion 
flow in both these states is uncertain, and there is 
no reason to assume complete alignment between the 
accretion flow angular momentum and black hole spin 
axes in these states. Given the time-dependent nature 
of the ray tracing, we can analyze the variability of the 
simulated accretion flows for the simplistic emission 
models used here to analyze the shape of their power 
spectra and to look for possible QPOs. 

The best time sampling of the simulations is in 90h and 
915h, which are used here at 8 observer azimuths, 3 incli- 
nations and 3 observed photon energies using the thermal 
emission model. Each light curve captures roughly 6 (20) 
orbits at r = 25M (10M), corresponding to a total ob- 
server time, Ai bs = 0.23Mio S, where Miq is the black 
hole mass in units of 10M Q . This duration is about 1/8 
of the total precession period for the torus in the 915h 
simulation. 

Figure [9] shows sample light curves and power spec- 
tra from the thermal emission model at lOkeV for an 
observer inclination of i = 60°. The secular trend is re- 
moved by subtracting the linear best fit from the light 
curve before computing the power spectrum. 

All power spectra are well fit by broken power law mod- 
els of the form: 



P{v)=Av~~ tl v<v b 



(1) 



where 71, 72 are power law indices and the break fre- 
quency u b , lies near 100Hz in both simulations. 
The tilted disk power spectra tend to flatten out at the 
highest sampled frequencies, ^1000Hz M^ . Figure fTOl 
shows median power spectra for the three different incli- 
nation angles from each simulation. The error bars are 
estimated from the standard deviation in log Power over 
observer azimuths and photon energies. At higher incli- 
nations, the peaks in the power around 100Hz grow, espe- 
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Time (s) Time (s) Frequency (Hz) 

Figure 9. Sample light curve and linear fit (left), light curve with linear fit subtracted (middle) and power spectrum (right). The units 
are scaled to a 10 solar mass black hole. 




Figure 10. Median power spectra for i = 30°, 60°, 90° from the 90h and 915h simulations. The errors are estimated from the standard 
deviations of the set of power spectra at observed photon energies of 1, 3, 10 keV at four observer azimuths. All power spectra are well 
described by the broken power law model, with break frequencies around 100Hz. 



cially for the tilted simulations. This would be expected 
from a source of excess power in the inner radii, where 
the larger Doppler shifts at higher inclination would en- 
hance the signal. 

To quantitatively compare the power spectra between 
the two simulations, the ratio between median power 
spectra in untilted and tilted simulations is plotted for 
each inclination in Figure 1111 The values are normalized 
to the combined uncertainties at each frequency. The 



overall plots arc shifted according to the mean ratio be- 
tween power spectra. 

At almost all frequencies, these ratios are within 
±2cr, and are unlikely to be observed as significant fea- 
tures. However, there are a few noteworthy features near 
lOOMfu H z. These are particula rly interesting given the 
finding by IHenisev et al.l (|2009T ) that the tilted simula- 
tion 915h contains excess power due to trapped inertial 
waves at mM^Hz. 
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Table 2 

Broken Power Law Fit Parameters 



100 

Frequency (Hz) 



1000 



Figure 11. Difference in the logarithms of median 90h and 915h 
power spectra, normalized to their combined standard deviations 
for i = 30°, 60°, 90°. The median 90h power spectra are shifted 
to account for their lower mean power. 




100 

Frequency (Hz) 

Figure 12. Sample power spectrum (solid), best fit broken power 
law model (dot-dashed) and upper and lower 99% (dotted) and 
99.9% (dashed) significance contours. 

To assess the significance of possible features in the 
PSDs, the power spectrum is fit with a broken power 
law model. The parameters from the best fit are used 
to simulate many random light curves with the same pa- 
rameters, and which contain no significant features. The 
significance is determined by comparing the values for 
the power at each frequency for each model power spec- 
trum with the distribution of random ones. An example 
is shown in Figure where a single power spectrum 
from the 915h simulation is shown, as well as the best 
fit broken power law model and upper and lower 99.9% 
confidence intervals from simulating random light curves. 

No obvious QPO features show up in this analysis. 
In several of the 9I5h light curves, the feature near 
SOM^Hz shows up as 99.9% significant. It appears at 
high significance in more of the light curves at high incli- 
nations. In the 90h simulations, almost all significant fea- 
tures are found at very high frequencies ~1000M-[^ Hz. 
These are spurious, caused by slight errors in the fit to 
the post-break slope incurred by ignoring all frequen- 
cies larger than 800i\/ 10 1 Hz. Including the highest fre- 
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Figure 13. Fractional difference in shell-averaged angular mo- 
mentum between tilted and untilted simulations with similar spins. 
The untilted simulations are nearly geodesic, while the tilted sim- 
ulations are increasingly sub-geodesic with decreasing radius. 

quencies in the fit can favor models with break frequen- 
cies ^500 Af 10 1 Hz, steep initial slopes and shallow post- 
break slopes. This occurs due to the denser sampling 
of the PSD at high frequencies. Simply ignoring the 
highest frequencies gives better results than a variety 
of more complicated weighting schemes. The features 
near 100Af{ 1 Hz from Figure [T^] never show up at more 
than 99% significance. In general, while the feature near 
SOAf^Hz in the tilted simulations is more convincing 
than anything from the untilted simulations, it does not 
appear at high enough significance at enough observer 
frequencies and azimuths to be identified as a QPO. 

Finally, fitting the sets of power spectra provides a 
general idea for the range of best fit values of the broken 
power law parameters. The median parameters found 
from the tilted and untilted simulation are listed in Ta- 
ble [2j where the quoted uncertainties are the standard 
deviations from light curves with different observer az- 
imuths and frequencies. Break frequencies have the units 
M^Hz. The break in slope becomes more pronounced 
at higher inclination as the initial slope becomes shal- 
lower while the post-break slope becomes steeper. The 
post-break slope is slightly shallower in the tilted simula- 
tions, while the initial slope is more strongly dependent 
on inclination in the untilted case. 

4. PHYSICAL CAUSE OF DISK TRUNCATION 

The observable signatures of tilted disks discussed so 
far are, for the most part, due to two main differences 
between tilted and untilted disks: tilted disks precess , 
and they are truncated outside r ms . iFragile et al.l (|2007l ) 
already discussed why the simulated accretion flows pre- 
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Figure 14. Shell-averaged entropy distributions for all simu- 
lations. Excess entropy inside r ~ 10M is generated by non- 
axisymmetric standing shocks in the tilted simulations. 

cess. It is our interest to better understand the physical 
cause for the large truncation radius. 

The first thing to note is that rapidly rotating black 
holes provide more centrifugal support to an accretion 
disk than slowly rotating black holes. Therefore, the an- 
gular momentum extraction mechanism at play in the 
tilted disks must be more effective at higher spin. This 
is confirmed in Figure 1131 where we plot the difference 
in density-weighted, shell-averaged specific angular mo- 
mentum for tilted and untilted simulations of comparable 
spin. The angular momentum is defined as t — —u^/ut, 
where is the fluid four-velocity and the shell-average 
of a quantity x is given by, 



(2) 



where f2 is the coordinate solid angle, g is the metric 
determinant and A — J J dQ^/—g. The density- weighted 
shell-average of x is defined as (pa;) /(p). The angular mo- 
mentum profiles for the untilted simulations are nearly 
geodesic outside of r ~ 5M. Inside of r ~ 10M, the tilted 
simulations become increasingly sub-geodesic, with the 
higher spin cases deviating more than the lower spin ones. 
The same trend holds when comparing the tilted simu- 
lations to the analytic result for the angular momentum 
profile of material on geodesic orbits in an equatorial disk 
i nclined 15° to th e black hole spin axis. 

iFragile fc Blaesl (|2008[ ) suggested that the non- 
axisymmetric standing shocks that occur in the inner 
radii above and below the midplane of the disk may 
enhance the outward transport of angular momentum, 
causing fluid to plunge from outside the marginally sta- 
ble orbit. To connect the enhanced angular momentum 
loss of the tilted disks with the standing shocks, we next 
look at a plot of the density- weighted, shell-averaged en- 
tropy profiles in Figure 1141 Since these simulations con- 
serve entropy except across shocks, the excess inside of 
r ~ 7M in the tilted simulations signifies the presence of 
extra shocks. The steepness of the entropy gradient gives 
some measure of the strength of these shocks. Again, we 
see that the effect is greatest in the simulations with the 
fastest spinning black holes. 

Further evidence linking the sub-geodesic angular mo- 



mentum profiles of the tilted simulations with the stand- 
ing shocks can be found from looking at the time- 
dependence of the shell-averaged angular momentum. 
While the untilted simulation remains nearly geodesic, 
the tilted simulations are continuously transporting an- 
gular momentum outward from r ~ 10M for the first 
~ 5000M before reaching a steady state, as would be 
expected from a dynamical mechanism. Finally, verti- 
cally integrated contour plots such as Figure [15] show that 
the angular momentum in the tilted simulations is non- 
axisymmetrically distributed. The regions of depleted 
angular momentum correspond to the standing shocks, 
which appear as regions of excess entropy in the bottom 

panels of Figure [15| 

Following IFragile fc Blaesl (|2008f) . we postulate that 
the standing shocks are caused by deviations from cir- 
cular orbits near the black hole. Figure [T5] shows the 
shell-averaged eccentricities of the orbits in each simula- 
tion, estimated at one scale-height in the disk using 



r 9(/3sin7) 



6M 



dr 



(3) 



where /3 is the tilt and 7 is the precession of each or- 
bital shellQ All quantities are calculated from fitting 
th e shell-averaged dis k tilt and twist (Eqs. 32 and 41 
of IFragile etafl 120071) with power laws, and using the 
resulting expressions in Equation ([3]). The increase in 
eccentricity toward smaller r adii leads to a crowding o f 
orbits near their apocenters (jlvanov 
which leads to the formation of the standing shocks. The 
eccentricity is larger for higher black hole spin, except in- 
side the plunging region where the fits become poor and 
the eccentricity is ill-defined. Equation ([3]) may indicate 
how these results depend on the initial tilt of the simu- 
lations. If we assume that the strongest dependence of e 
on tilt is through /3 and that d/3/dr and dj/dr remain 
unchanged for different tilts, then equation (|3]) suggests 
that the eccentricity of the orbits should vary roughly 
linearly with the initial tilt, at least for small angles. 
This prediction is tentatively confirmed by a simulation 
we have done that started with an initial tilt of 10°. 

5. DISCUSSION 

Tilted accretion flows will inevitably be present in a 
significant fraction of black hole sources with L/L ec id J$ 
0.05 and possibly L/L e dd ^ 0.3 (thick or slim disks). Us- 
ing relativistic ray tracing and a set of simple emissivities, 
we have compared the radiation edge, emission line pro- 
files and power spectra of simulated black hole accretion 
flows with a tilt of 15° to their untilted counterparts. We 
find the radiation edge is independent of black hole spin, 
while the untilted simulations agreed with the expected 
qualitative trend of decreasing inner radius with increas- 
ing spin . These result s for the radiation edge confirm the 
work of iFragili (|2009f ) . who used dynamical measures to 
locate the inner edge. 

Due to the independence of inner edge on spin, the 
red wing of tilted accretion flow emission line profiles is 
also fairly independent of spin. This introduces a possi- 
ble complication for attempts to measure black hole spin 



1 This definition of e differs from Ivanov & Illarionov ( 1997) by 

Blaesi H2008I 1 used the formula 



a ph ase factor of 7r/2 in 7 
from [Ivanov & Il larionovl 
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1997) without modification. 
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Figure 15. Vertically integrated contour plots of specific angular 
momentum (top) and entropy (bottom) for snapshots of the 90h 
(left) and 915h (right) simulations. The color scale is linear, in- 
creasing from blue to red to yellow to white. The non-axisymmetric 
shocks in 915h correspond to regions with deficit (excess) angular 
momentum (entropy). 
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Figure 16. Shell-averaged orbital eccentricities for all tilted sim- 
ulations, estimated at one scale-height in the disk. The increasing 
eccentricity of orbits toward smaller radii leads to a crowding of 
orbits at their apocenters, which, in turn, can generate standing 
shocks. 

from sources which may be geometrically thick. In gen- 
eral, measurements of small spin (large inner radius) may 
be unreliable unless the disk is known to be untilted. A 
reliable estimate of a large black hole spin (small inner 
radius), in contrast, could rule out the presence of a tilted 
disk. 

The blue wing can be much broader for tilted ac- 
cretion flows, and the tilted-disk line profiles depend 
strongly on the observer azimuth as well a s inclination. 
Since a tilted disk is expected to precess ([Fragile et al.l 
I2007T ) . highly variable emission line profiles could signify 
the presence of a tilted a ccreti on flow, as pointed out 
by lHartnoll fc Blackman (2000) for warped thin disks. 
Since many LLAGN and X-ray binaries in the low/hard 
state should be tilted, time- variable emission lines should 



be quite common, and this effect is unlikely to signifi- 
cantly depend on accurate reflection spectrum modeling. 
Although the simulations can only be run for a short 
time compared to the precession time scale, precession 
is a possible source of low frequency quasi-periodic oscil- 
lations when the accretion flow is optically thin due to 
the modulation of Doppler shifts as the velocities in the 
accretion flo w align and misalign with the observer's line 
of sight (see llngram et al.l [2009. for more discussion of 
QPOs from precessing tilted disks). 

Finally, we have studied power spectra for our 
simple models. We find broken power law spec- 
tra with break frequencies around 100A/ 10 1 Hz and 
power law indices in the range 0-2 (3-4) pre- 
(post-) break for both tilted and untilted s i mula- 
tions. Previous stu dies (fArmitag e fc Reynolds! 120031 : 
INoble fc KroTikl [20091) found si ngle p ower laws with in- 
dex ~2. lArmitage fc Reynolds! (|2003l) found that power 
spectra from individual annuli are well described by bro- 
ken power laws where the break frequency is close to 
the local orbital frequency - the averaging of many an- 
nuli with an emissivity that falls with radius smooths 
the power spectrum into a single power law. We see 
the same behavior in our simulations; the break frequen- 
cies from power spectra of individual radial shells agree 
with the local orbital frequency for both simulations 90h 
and 915h. A break frequency lOOAT^Hz then implies 
a radius of r 16M. Our broken power law spectra are 
therefore likely due to the fact that our emissivity peaks 
relatively near the outer radius used for the ray tracing, 
r = 25M. A larger radial domain would likely shift the 
break to smaller frequencies. 

Observed break frequencies in the low/hard state are 
typically i/f, ~ 0.1 — 1Hz, which may be caused by the 
transition from a thin disk to a thicker, ADAF flow 
(E sin et al.lll997l ). That would imply a transition radius 

r t ~ 200 - lOOOM-^M. Our results for pre- and post- 
break slopes from both tilted and un tilted simulations 
agree with those found in Cygnus X-l ([Revnivtsev et al.l 
[2000h for an inclinati on i = 30°. In GRO J1655-40 
([Remillard "et~aT1[l99"9l) our pre-break slopes agree for all 
inclinations. However, the PSD for that source is well 
described by a single power law. 

There is no clear evidence in our work for high fre- 
quency QPOs due to t he tra pped inertial waves identi- 
fied by iHenisev et al.l (|2009f ). although there are more 
features in power spectra from the 915h simulation at 
higher significance than in 90h. Even when computing 
PSDs for sets of spherical shells from the simulations, 
there are no clear features in the tilted power spectra 
that are not also present in the untilted case. It is possi- 
ble that this result could depend on the chosen emissiv- 
ity. Alternatively, the excess power in trapped inertial 
waves could be insufficient to rise above the red noise 
continuum. 

The independence of the inner radius of the tilted sim- 
ulations on black hole spin is attributable to the extra 
angular momentum transport provided by the asymmet- 
ric standing shocks. These shocks are only present in the 
tilted simulations. Their strength scales with black hole 
spin, which is a necessary condition for countering the 
greater centrifugal support at higher spins. The stand- 
ing shocks, in turn, appear to be attributable to epicyclic 
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motion within the disk driven by pressure gradients as- 
sociated with the warped structure. Again, this effect 
scales with the spin of the black hole, which contributes 
to the stronger shocks. 

For small tilt angles, the orbital eccentricity scales as 
e *~ /3. This suggests that significant deviations be- 
tween the spin-dependence of the radiation edge and the 
marginally stable orbit should be present even at mod- 
est tilt angles j3 > 5°. At larger tilts, it is unclear if the 
increasing eccentricity will lead to an inner edge that in- 
creases with spin. This is both due to the uncertainty in 
the radial tilt and twist profiles j3(r) and j(r) at larger 
tilts, and to the lack of a quantitative connection be- 
tween inner di sk edge and ec centricity. The dynamical 
measures from iFragild ([2009D place the location of the 
inner edge in a simulation with a = 0.9M and j3 = 10° 
closer to the location of 915h than 90h. This data point 
supports the idea that a noticeable departure between 
r c dg C and r ms should exist between tilted and untiltcd 
disks even for (3 > 5°. It also suggests that at larger 
tilt angles, r c d gc is likely to increase with spin unless the 
effect saturates at /3 rj 15°. Simulations with larger tilt 
angles will be able to address this question with certainty. 
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